Spin-orbit coupling in interacting quasi-one-dimensional electron systems 
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We present a new model for the study of spin-orbit coupling in interacting quasi-one-dimensional 
systems and solve it exactly to find the spectral properties of such systems. We show that the 
combination of spin-orbit coupling and electron-electron interactions results in: the replacement of 
separate spin and charge excitations with two new kinds of bosonic mixed-spin-charge excitation, 
and a characteristic modification of the spectral function and single-particle density of states. Our 
results show how manipulation of the spin-orbit coupling, with external electric fields, can be used 
for the experimental determination of microscopic interaction parameters in quantum wires. 
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In contemporary condensed matter physics there exists 
a great variety of electron systems that can be considered 
as quasi-one-dimensional (QID). Among them are con- 
ducting polymers |l| , carbon nanotubes |^] , and semicon- 
ductor heterostructures [^. They combine the richness 
of observable physical properties with the possibility of 
finding exact solutions to non-trivial interacting prob- 
lems. The major theoretical advance in this field was 
the formulation and solution of the Tomonaga-Luttinger 
model which revealed features that are generic to 
many interacting QID electron systems, such as sepa- 
ration of charge and spin degrees of freedom and the 
anomalous scaling of correlation functions |^ . 

The very existence of QID systems (or quantum wires) 
arises essentially from confining electrical forces, which 
prevent the lateral motion of electrons. The accompa- 
nying electric fields necessarily give rise to a coupling 
between the spin and orbital degrees of freedom of prop- 
agating electrons — the spin-orbit (SO) interaction. Al- 
though this interaction is relativistic in nature, it never- 
theless results in a significant modification to the band 
structure of QID systems (see, e.g., experimental |^ and 
theoretical [Q papers). In strictly ID and 2D systems, 
the role of the SO interaction basically reduces to a "hor- 
izontal" (in momentum space) splitting of the energy 
branches corresponding to spin-up and spin-down elec- 
tron states. In contrast to this, the existence of quan- 
tised transverse energy subbands in quantum wires brings 
about an additional important manifestation of the SO 
coupling: the absence of a vertical symmetry axis for 
each spin branch (see Fig. ^ and Ref. [|j ) . This absence 
results in a breakdown of chiral symmetry so that the 
electron Fermi velocities become dependent on the direc- 
tion of motion. This effect monotonically increases as 
the SO coupling is enhanced. As applied to interacting 
systems, the absence of chiral symmetry raises two fun- 
damental questions: how does the system respond to the 
asymmetry of the single-particle spectrum; and what is 
the fate of the spin-charge separation in the presence of 



the SO coupling? It is the aim of this Letter, to establish 
a realistic model for the SO coupling in interacting ID 
systems which answers these questions, and to provide 
experiment with a set of quantitative predictions which 
will allow their verification. 
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FIG. 1. Single-electron energy spectrum in a quantum wire 
with spin-orbit interaction. The dotted line is the Fermi en- 
ergy level E = Ef and ±pi,2 are the Fermi momenta of re- 
spective groups of electrons. 



In constructing a Hamiltonian for such a system, we 
consider the case where the Fermi energy Ep is suffi- 
ciently small, such that only the lowest energy subband 
in a quantum wire is partly filled, while all the others are 
empty. This regime has proven to be the richest in non- 
trivial experimental results [^,^, and the SO effects are 
expected to be most pronounced here. As a natural 
way of capturing the essential physics in a quantum wire, 
we suggest the use of a modification to the Tomonaga- 
Luttinger model that takes into account the asymmetric 
single-particle spectrum in Fig. |^. Namely, we take the 
Hamiltonian to have the form H = Hq + Hint , where Hq 
describes the linearized spectrum, 
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-iv2 J dx [Tp^^ ^^d^ipR^i - V'IjC'^'/'l.t) , 
and Hint is responsible for the EE interaction: 
Htnt = 92 j dx '0lj,tV'fl,TV'l,xV'L,i 

+ 94, dx '0fl,t'0i?,TV'lj,4V'i?,i + (i? ^ . 



(1) 



where 



(2) 



The operators ^r,s{x) {r = R, L; s =t, i) annihilate spin- 
up (t) and spin-down (|) electrons near the right (i?) and 
left (L) Fermi points. In Eq. (Q) the electron velocities 
vi,2 = dE{pi^2) / dp are derived from the dispersion law 
in Fig. 1^. In general, vi ^ V2 as long as the spin-orbit 
interaction is finite 

In the interaction Hamiltonian (H) we include only for- 
ward scattering contributions, while the backward and 
Umklapp scattering are left out. Neglecting the Umklapp 
scattering is known ||] to be safe if energy bands are far 
from being half- filled, which is exactly the case, e.g., in 
quantum wires patterned in semiconductor heterostruc- 
tures ||. As regards the backscattering, it can be shown 
that, similar to the case of zero SO coupling it renor- 
malizes to zero for repulsive interactions and therefore 
can also be omitted. In this letter, we concentrate on the 
major tendencies arising from the SO coupling in inter- 
acting electron systems, therefore we assume momentum- 
independent interaction potentials 52 > and (74 > in 
Eq. (H), which corresponds to a well-screened repulsion 
between electrons. 

The Hamiltonian (|l|), (H) is reminiscent of that for 
the multi-component Tomonaga-Luttinger model [^js) 
which consists of mutually interacting Luttinger liquids 
with different Fermi velocities. However, in contrast to 
our case, the model [|o[||] assumes that each liquid has 
a symmetric single-electron spectrum (as occurs, for ex- 
ample, with Zeeman splitting [pT[| ) and therefore, as we 
will show, it describes qualitatively different physical be- 
haviour. 

We study the interacting model (||) with the help 
of the bosonization technique Within this for- 

malism, fermionic operators are represented as 

'^r,s{x) ^ exp(— ■i^'rjs), where the operators '^r,s{x) are 
linear combinations of so-called phase fields ^p{a) {x) and 
^p{a){x) that obey the bosonic canonical commutation 
relations (see Refs. |12 for details). In terms of these 
phase fields, the Hamiltonian (^, (||) takes the bosonized 
form: 

H^^jdx[ VpKp{T:npf + ^^{d.Vp? 

+ Sv [(7rHp)a,vJ, + (7rH„)<9,(^p]}, (3) 



"0 = (vi + V2)/2, Sv = V2- Vi 



1 /2 

and Vp(^„) = [{vo ± gi/^irf - {g2/2-Kf] and Kp(„) = 

1 /2 

[(27rz;o T 52 ± 54)/(27rz;o ± ^2 ± 34)] ■ The first (sec- 
ond) line in Eq. (^) describes free propagation of charge 
(spin) density wave with the velocity i'p((t) and "stiff- 
ness" coefficient Kp{^„y The third line in Eq. (^) is 
proportional to the velocity difference 5v and therefore 
represents the strength of the SO interaction. In the 
Tomonaga-Luttinger model, where vi = V2, this term 
is absent and this results in a decoupling of the charge 
and spin degrees of freedom: the so-called spin- charge 
separation [Q. This decoupling inhibits the existence of 
quasi-particles with spin 1/2 and charge — e, the basic 
excitations of a Fermi liquid, and gives rise to a different 
state of matter, the Luttinger liquid, which has bosonic 
excitations in the form of independent spin and charge 
density waves. 

As the SO interaction is switched on, i.e. at 6v ^ 0, 
the third term in Eq. (||) starts to affect the dynamics 
of the system. It couples together the p- and cr-fields 
and thereby destroys this spin-charge separation. A sim- 
ilar effect is also found in a spinful Luttinger model in 
a magnetic field In our case, the asymmetry of 

the single-electron spectrum in Fig. |l| results in a differ- 
ent mechanism for the violation of spin-charge separation 
from Ref. We stress however, that our model is ex- 
actly solvable since the Hamiltonian ([l]), (§) is quadratic 
in the phase fields. Here we choose to solve it within the 
functional integration formalism in imaginary (Matsub- 
ara) time (Tsj. 

The most profound effects of the SO coupling on in- 
teracting QID systems can be seen in the behaviour of 
single-particle characteristics, such as the spectral func- 
tion Pr{q, i^), 
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r = i?, L, 



(5) 



and the density of states TV (w) = (l/27r) / ^1 Pr{q^^)- 
Here q and u> are the momentum and energy of elemen- 
tary excitations and Gi^J*\q,uj) is the Fourier transform 
of the retarded Green function 



Girf{x,t) = ie{t) ({^,,,(x,t), vJ,.(o,o)}) 

= -ie{t) [Gr,s{x, t) + Gr,si-X, -t)] 



(6) 



where 9{t) is the Heaviside step function and Gr.s{x, t) — 
(Vr,s(a;, T)'(/'r ,,(0, 0)) 1^^.^ with r being the Matsubara 
time. Leaving detailed calculations to a regular article, 
we give here an exact expression for Gr,s{x,t) at T = 0: 
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Gr,s {x, t) 



A 



A + i(uit — x) 



A 



A + i{uit + x) 



A 



A + i{u2t — x) 



A 



(7) 



A small (of the order of the lattice parameter) length 
scale A defines ultraviolet cut-off From Eq. (|^) it 

follows that although the SO coupling destroys the spin- 
charge separation, it nevertheless preserves the anoma- 
lous scaling of correlation functions. At Sv = 0, Eq. (|^) 
coincides with the well-known expression ||^. The expo- 
nents of {i = 1,2) are determined by values of 'Up((T), 
Kp[a), Sv and depend on the chiral r and spin s indices. 
For arbitrary interaction constants 52 and 54, the expres- 
sions for of are rather lengthy and here we give them for 
the realistic case of equal interaction strength between 
all the spectral branches, i.e. for 92 = .94 = {t^vq)P: 

0fir,s) = {-iycfir,s)r^,/2irjl~-rjl), (8) 

C±(r, s) = ±r [(1 + rse/2)V77, - r,,] - (1 - rse/2) 

+ [(l + r.6/2)(l-6V4)-/3V2] hi 

Here e = 6v/vq, and on the r.h.s. of the last equation 
r = ±1 for right (left) spectrum branches and s = ±1 for 
spin-up (spin-down) electrons. 

The quantities 771^2 = '7i,2(e) defined by 



(9) 



are dimensionless velocities of the new independent 
bosonic excitations that replace the spin and charge den- 
sity waves, respectively, as the single-electron spectrum 
becomes asymmetric (e ^ 0). Each of these new ex- 
citations are superpositions of the old ones, they carry 
both charge and spin, and have a sound-like spectrum 

= ui^q. For e = 0, we have wi,2 = Wo-(p) and return 
to the spin-charge separation. Eq. (H) demonstrates that 
increasing e at constant (3 pushes ui and U2 away from 
Va and Vp as well as from each other: one of the exci- 
tations monotonically accelerates (u2 grows) while the 
other monotonically slows down (ui decreases). This ef- 
fect becomes more pronounced in systems with stronger 
EE interaction. Since the parameter e is controlled by 
the SO interaction and therefore by the lateral electri- 
cal confinement, the strength of this confinement must 
affect drastically the dynamics of elementary excitations 
in quantum wires. Further increasing e eventually re- 
sults in a vanishing of the velocity ui at the "critical" 
point eo — 2^/1 — (3 and in a possible phase transition 
accompanied by phase separation This phenomenon 
is beyond the scope of this Letter. 

Eq. (0) allows one to calculate the spectral function 
Pr{q,uj) which can be obtained in magnetotunneling ex- 
periments [ p^yr^ ]. Fig. H shows the spectral function 
pii(q,u}) for right-moving electrons in the range of en- 
ergies Lu where all singularities of ppi{q,uj) are located. 



At each value of e < eo, the spectral function has two 
singular points corresponding to the bosonic excitations 
with velocities ui and U2 [see Eq. (^]. In full accor- 
dance with Eq. (^, the distance between the singular 
points grows with e. The singularities of pfl(g, w) are 
governed by power laws with exponents depending on 
Wp(o-), -K'p((t)7 and e. The dependence on e is such that 
one of the singularities becomes sharper while the other 
becomes smoother as e grows. It is noteworthy that the 
dependence of or on q at fixed to is very similar to that 
shown in Fig. g and exhibits the same tendencies as a 
function of e. 
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FIG. 2. Spectral function pR{q,uj) (in arbitrary units) for 
a fixed q = O.OIA"^ and several values of e and for /3 = 0.5 
Wo /A). Note that pL{q,w) = pR{q, -w). 
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From Eqs. (||), (||), and (|^) one can also calculate the 
single-particle density of states N{u}). This quantity may 
be obtained in tunneling measurements since the tunnel- 
ing current, e.g., between a wide metal and an interacting 
quantum wire is proportional to N{eV), where V is an 
applied voltage (see Ref. |l^]). It can also be determined 
in angle- integrated photocmission experiments [p^ . For 
uj —^ we obtain the following behaviour: 



N{uj) 

Nn 



1 



{uj/uja) 



i+e2-i 



if^-rf- r(0i+02) 



= &f-rO-. (10) 



Here = v^j A is the natural energy imit of the order 
of the Fermi energy, iVo = 2 /ttuq is the constant density 
of states at /3 = e = 0, and V{x) is the Gamma function. 

Fig. ||a demonstrates the effect of the EE interaction 
on the density of states N[uj) at a fixed strength of the 
SO coupling. At zero EE interaction (/? = 0), the func- 
tion is a constant, iV(w)/iVo = 1/(1 - e^/4). As 
the EE interaction is switched on (/3 7^ 0), it starts to 
affect the formation of the lowest lying excitations and 
N(lo) takes on a power-law behaviour in the vicinity of 
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u) = 0. The width of this power-law interval becomes 
progressively larger as (3 grows and leads to a monotonic 
suppression of the density of states. This effect is also 
present for zero SO coupling ||^. 
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FIG. 3. Normalised single-electron density of states: (a) 
e = 0.5 and (3 = 0.0, 0.1, . . . , 0.9; (b) f3 = 0.6 and 
e — 0.0, 0.1, . . . , 1.1. The dashed curves correspond to e = 0.5, 
l3 = 0.6. 

Fig. ^ shows the evolution of the function N{uj) for 
fixed (3 as the SO coupling is varied. For large values of 
uj, where the role of EE interactions is not significant, 
the magnitude of 7V(a;) increases as e grows. However, 
at small values of uj, where the nature of the elementary 
excitations is essentially dictated by EE interactions, the 
effect of the SO coupling on N{u!) is qualitatively the 
same as that of the EE interaction, that is, increasing e 
leads to a suppression of the density of states. We empha- 
size that in the standard multi-component Tomonaga- 
Luttinger mod el ||l0| ,|5|] as well as in spin-polarized Lut- 
tinger liquids there is no interval of u) where the 
density of states is suppressed by increasing the veloc- 
ity difference between spectral branches. The existence 
of such an interval is a unique manifestation of the SO 
coupling in quantum wires. 

The characteristic dependence of the spectral function 
p{q,Lu) and the density of states N{uj) on the microscopic 
parameter e, for fixed /?, should enable one to extract 
both quantities from experiment and thus determine how 



strong the EE and SO interactions are. Such dependen- 
cies can be obtained experimentally by changing the SO 
coupling directly by varying an electric field perpendic- 
ular to a quantum wire (quantum-well field) as it was 
done, e.g., in Ref. [Q. 

Theoretical calculations Q of the electron bandstruc- 
ture modified by the SO coupling indicate that the value 
of e in typical QID semiconductor systems should be 
~ 0.1 - 0.2 and appears sufhciently large to observe the 
principal tendencies in the behaviour of p(g, lo) and N{uj) 
caused by the SO couphng. 

In conclusion, we have formulated and solved analyt- 
ically the problem of the interplay between EE and SO 
interactions in QID electron systems and found that: (i) 
the spin-charge separation breaks down while the anoma- 
lous scaling of correlation functions is preserved; (ii) the 
single-particle characteristics, such as the spectral func- 
tion and the density of states, are essentially modified 
and controlled by the strength of the SO coupling; (iii) 
varying the SO coupling with the external electric field 
can be used to extract the microscopic parameters of 
quantum wires. 
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